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VIBRATIONS OF THE UNIVERSE 

by 

E. R. Harrison 

ABSTRACT -Q f) $ Q 

The aim is to review and extend our knowledge of the time de- 
pendence of fluctuations in density in homogeneous and isotropic 
models of the universe. Perturbations are assumed to be small and 
a linearized, normal mode analysis is used. Only gravitational in- 
teractions and irrotational motions are considered. Using first 
Newtonian and then general relativity theory it is shown that the 
amplitude of various modes increases in expanding and contracting 
models of the universe. The origin and growth of celestial struc- 
ture require that uniform models of the universe are unstable 
against arbitrarily small perturbations. However, the rate of growth 
of the various modes, particularly in an expanding universe, is suf- 
ficiently slow to cast doubt on the instability of the models. A lin- 
earized gravitational theory appears inadequate to account for the 
origin of structure in the universe; hence alternative explanations 
must be found. f\ -■*- / 
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VIBRATIONS OF THE UNIVERSE 


*■ 


1. INTRODUCTION 

Nobody believes that the universe was created in its present form, complete 
in all its array of celestial detail. If the inchoate universe is featureless, then 
we must show how differentiation and structure contrive to evolve. The physics 
and evolution of stellar and galactic structure is the concern of astrophysics. 

But the nature and development of an environment favorable to the formation of 
astrophysical objects is the concern of cosmology. 

The study of the vibrations of the universe is a fascinating subject, apart 
from any other consideration. However, in an attempt to account for the actual 
universe, as distinct from the smooth idealized models of cosmology, we might 
hopefully expect that the universe is unstable. In this way, various modes of 
vibration of arbitrarily small amplitude grow relatively rapidly and eventually, 
in a time short compared with the age of the universe, the foundations of astro- 
physical structure are laid down. However, it turns out that in an expanding 
universe the growth rates of the various modes are too small to account for 
appreciable irregularity. 

We assume that all modes have small amplitudes and use linearized equa- 
tions. Only gravitational interactions are considered; the motions are assumed 
to be irrotational; and a simple equation of state is used. The behaviour of the 
modes is then studied, first with Newtonian and then general relativity theory. 
Some general comments are made in Section 6. 
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2. ORIGIN OF STRUCTURE IN THE UNIVERSE 

2.1 Initial Conditions 

Cosmology is the science of the universe; at present it seeks only to account 
for the macroscopic nature of the universe. Instead of the physical universe, 
cosmology deals with a featureless, idealized universe that is everywhere iso- 
tropic and homogeneous and contains a uniform fluid of rudimentary properties. 1,2 
Out of such generalizations emerge several simple models, and one's choice is 
largely an act of faith. Cosmology, in spite of its antiquity, is in an immature 
state; ultimately it must bridge the gap between its present idealizations and the 
physical universe with realistic models. 

The provision of background and initial conditions for the origin and forma- 
tion of structure is a cosmological problem. The subsequent evolution of struc- 
ture into its detailed manifestations lies in the provinces of cosmogony, astro- 
physics, and every other science. The physical universe, as distinct from the 
current cosmological models, is complex and diverse. From the simplest of all 
possible points of view the diversity consists of variations in the density and 
motion of matter. Given a cosmological model that is a valid description of the 
universe in the large, it should be possible to show, with refinements of the 
cosmic fluid when necessary, that perturbations are capable of evolving in time 
into configurations of density and motion which resemble the grosser features 
of the physical world. The first step therefore is to inject some realism into the 
cosmological models by perturbing their fluid density and motion. The behaviour 
of the perturbations might then provide a means of judging whether a model is 
acceptable or not. 
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Two hypotheses are possible concerning the initial conditions in an evolving 
universe:* 

A. Structure is created with the universe, or imprinted in it from the earliest 
moment of its expansion. 

B. Structure develops in the universe from initially small random 
disturbances. 

Hypothesis A is as old as cosmology. All current models (in which the cosmo- 
logical term is zero) evolve either from or through a singular state. The uni- 
verse originates before, at the time of, or after the singular state. It is impos- 
sible to imagine any structure surviving the singular state and therefore hypothesis 
A demands that the origin occurs after the singular state. Alternatively, by modi- 
fying the laws of physics , 4 the singular state can be avoided and the universe 
contracts to and expands from a finite size. This theory has the advantage that 
some structure might survive passage through the 'bounce' and act as dense 
nuclei around which further condensations subsequently form. It has the further 
advantage that in the linear theory, as we shall show, the growth rate of disturb- 
ances is usually larger in a contracting universe than in an expanding universe. 

The concept of condensations growing in a uniformly dispersed medium 
(hypothesis B) is as old as the idea of gravitation. Jeans points out that in some 
correspondence Newton remarks : 5 'But if the matter were evenly disposed 
throughout an infinite space, it could never convene into one mass; but some of 


These hypotheses do not cover the cose of the steady state model which is outside the range of 
this discussion. 5 
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it would convene into one mass and some into another, so as to make an infinite 
number of great masses, scattered great distances from one to another through- 
out all that infinite space. And thus might the sun and fixed stars be formed, 
supposing that matter were of a lucid nature.' In his own work on gravitational 
instability Jeans writes : 6 'We have found that, as Newton first conjectured, a 
chaotic mass of gas of approximately uniform density and of very great extent 
would be dynamically unstable; nuclei would tend to form in it, around which the 
whole of the matter would ultimately condense.' These comments were made 
with the idea in mind of a static universe. If we consider a universe already 
fragmented into widely separated 'islands,' each having a density large compared 
with the average density, such that the contents of each island do not partake in 
the expansion between the islands, then Jeans' comments are acceptable. But 
these are the initial conditions of cosmogony; the problem for cosmology is to 
explain how it is possible in the first place for islands to form. 

A principal difference between the two hypotheses is that in B perturbations 
of the universe can, at least initially, be studied with a linear theory, whereas 
in A it is doubtful whether a linear theory is valid at any stage in an expanding 
universe. 

In this discussion we shall discount ad hoc modifications of gravitational 
theory and hence A reduces to the proposition that the universe began with much 
of its detailed design predetermined. The obvious comment is that such ideas 
run counter to the main trend of science. By developing and employing the laws 
of physics, with their properties of symmetry and invariance, the aim is to 
rationalize our observations in terms of given initial conditions . 7 With the 
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advance of science the laws of physics account for more and more of the regu- 
larity of nature at the expense of the need to assume regularity in the initial 
conditions. One might hope that ultimately the laws of physics, operating in a 
universe containing random and irregular fluctuations of small amplitude, will 
be adequate to account for the observed macroscopic structure of the universe. 

The choice between hypotheses A and B is a matter of preference; if A, then 
there is scarcely any problem for cosmology; but if B, then it is necessary to 
show at least that in a linear theory certain fluctuations can grow in an expand- 
ing universe. In the following we adopt hypothesis B. 

2.2 Linear Stability Theory 

It is assumed that all disturbances are small and consist of a superposition 
of normal modes of a complete set; only gravitational interactions are considered 
using a linearized theory. 

Within a comoving system of coordinates the uniform fluid of the cosmologi- 
cal models is in a hydrodynamic stationary state. A system in a stationary state 
is unstable when a small disturbance grows in time and leads to a changed con- 
figuration of the system. Thus, if one or more modes is time-growing, and the 
characteristic growth time is short compared with the lifetime of the system, the 
system is unstable. 8 According to hypothesis B the universe is an unstable 
system. 

Planets, stars, stellar associations and clusters, galaxies, galactic clusters, 

form a hierarchy of perturbations in which the amplitude diminishes as its 

spatial extent increases. It follows that a linearized treatment is most valid for 
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fluctuations extending over extreme distances; that is for the lower vibration 
modes of the universe. As we go back in time the hierarchy of celestial objects 
dissolves into the increasing mean density of the universe, and the amplitude of 
all perturbations relative to the mean density diminishes. At a sufficiently early 
epoch, according to hypothesis B, the contrast density (ratio of the perturbation 
in density and unperturbed density) is small enough for a linear treatment. It 
follows therefore that a linear gravitational theory is cosmological in the sense 
that it is limited to small contrast densities, either remote in time, or extend- 
ing over cosmic distances. 

If the universe is unstable only for long wavelength perturbations we might 
conjecture with Jeans that there is a process of fragmentation 9 'of nebulae out 

of chaos, of stars out of nebulae, ,' and so on. Protogalaxies, or larger 

masses, first form because the universe is unstable. Regions of enhanced density, 
in which matter no longer expands with the universe, provide an environment 
more favourable for the formation of smaller condensations. Inhomogeneity, 
anisotropy, and complex properties of the fluid develop. No longer limited to a 
linearized treatment, we are free to invoke all the cosmogonic paraphenalia of 
turbulence, magnetic fields, radiation, dust, and so forth . 10 

Alternatively, if the universe is unstable for short wavelengths we might 
conjecture with Layzer that there is a process of clustering , 11 whereby small 
condensations first form and subsequently interact to create larger and larger 
self-gravitating systems. The arguments in favour of this process, however, 
have not progressed very far nor have they gained wide acceptance. 
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From the cosmological point of view the concepts of fragmentation and clus- 
tering are by no means mutually exclusive. The universe could be unstable for 
a large class of modes, or a wide spectrum of wavelengths, and the rate of growth 
of the different wavelengths determines whether elementary structure evolves 
by fragmentation or clustering, or by both processes acting simultaneously. Con- 
ceivably, early population II stars evolve out of inhomogeneities laid down either 
prior to or at the same time as those leading to galactic structure. All this, 
however, is speculation and we must wait for cosmology to give a clear and re- 
liable account of the origin of structure in the universe. 

3. NEWTONIAN COSMOLOGY 

3.1 Newtonian Models 

In 1934 McCrea and Milne 12, 13 used Newtonian theory to derive the equations 
of a universe obeying the cosmological principal. 2 It is assumed, as in general 
relativity, that in the unperturbed state there is everywhere a perfect fluid of 
uniform mass density p and isotropic pressure p. The equations are identical 
with those derived using general relativity theory, provided the pressure is negli- 
gible in comparison with the energy density pc 2 (c is the speed of light). Several 
authors 14 ' 18 have discussed the validity and limitations of Newtonian cosmology. 

As we shall show, the equations of Newtonian and general relativity theory 
for the perturbed state are similar when the pressure is negligible. The ad- 
vantage of the Newtonian treatment is its simplicity; furthermore, it provides 
physical insight which helps to reduce the general relativity equations to their 
simplest form. Before proceeding to the Newtonian equations of a universe in a 
perturbed state, the treatment for the unperturbed state is presented briefly. 
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Let r be the position at time t 0 of any element of fluid. At time t, let the 
position of the same fluid element be (S/S 0 ) r, where S(t) is a universal function 
of time and S Q = S(t Q ). This condition ensures that the fluid density p re- 
mains uniform and is a function only of time. Thus r is a comoving position 
vector and r , 6, <p are comoving spherical coordinates. The velocity and accel- 
eration of a fluid element are 


u 




( 1 ) 


where dots denote time differentiation. Within the comoving coordinate system 
the ordinary gradient operator V' transforms to (S Q /S)V. 

The equations of motion, continuity, and Poissions equation, are 


£ — = - V 4>* - I Vp , (2) 

S 0 dt P 


s 

s 


0 


dp 

V • (p u) , 

dt 


(3) 


(S Q /S) 2 V 2 i/>* = 4rrGp - A, 


(4) 


where p is the pressure, G the gravitational constant, and the gravitational 
potential. The cosmological term A is included for historical interest. In a 
uniform universe the Vp term in (2) vanishes. Using (1), and taking the diver- 
gence of (2), it follows that 
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3 S + (4-ttG/d - A)S = 0. (5) 

From the equation of continuity (3), 

pS 3 = constant, (6) 

and hence (5) can be integrated and becomes 

S 2 = I (8?7Gp + A)S 2 -k (7) 

3 

A universal constant of integration C is absorbed by the transformation: S - S | C | 
ands has now the dimension of time and k has the value of 1 , 0, or - 1 . 


Equations (6) and (7) are the Newtonian equations of an isotropic and homo- 
geneous universe. Furthermore, if 

S - c* 1 R, 

the equations are 


R 2 = I (8?tGp - A)R 2 - kc 1 
3 


( 8 ) 


pR 3 = constant, 


O) 


which are identical with those usually derived with general relativity, when the 
pressure is small compared with the energy density pc 2 (First derived by 


1/2 
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4 

Friedmann for p = 0, and 19 k = 1 and 20 k = -1.) The velocity of light is an alien 
intruder in Newtonian hydrodynamics and therefore we keep to (6) and (7) for the 
present. 

Using the constant 


fi = 8ttG / oS 3 /3, (10) 

the solutions of (6) and (7) are 


* 

II 

o 

S 

= /^y 2 ’ 


k - 1: 

S 

= sin 2 y , 

(ii) 

#c = - 1: 

s 

= (i sinh 2 y, 



and t(y) is found by integrating 


dt = 2 S dy . 

Within a sphere of radius r the potential energy is 


n = - 


16 

15 



7T 2 G p 2 r s 


/3 E 
kS ’ 


( 12 ) 


(13) 


and the kinetic energy is 


T = 


2 

5 


S 2 S 3 


s 5 

a o 


77 p r 5 


S 2 E 

K 


(14) 
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where E = Q + T is the total energy and A = 0. From (14): E = -*T/S 2 and the 
total energy is zero when k = 0. 

The advantage of the Newtonian equations is the ease with which they can be 
physically interpreted. Thus, assuming A = 0, when k = 0, E = 0, the fluid ele- 
ments have velocities equal to their escape velocity and their trajectories are 
of the parabolic class; and for k = 1 (k = -1), E < 0 (E > 0), the fluid elements 
have velocities less (more) than their escape velocity and their trajectories are 
of the elliptical (hyperbolic) class. These interpretations are less obvious when 
(8) and (9) are derived from general relativity. In that case * is the curvature 
constant and space is flat (k = 0), spherical or elliptical ( k = 1), and hyperbolic 
( * = - 1 ). 


3.2 Perturbed Newtonian Models 

Various authors 21 ' 29 have used Newtonian gravitational theory to study the 
time dependence of density fluctuations in a uniform fluid of finite or infinite 
extension. From the cosmological point of view a Newtonian treatment is scarcely 
adequate; not only is it limited to low density, but also the long wavelength modes 
of Euclidean space are inapplicable in curved space. Fluctuations at high density 
and large-scale fluctuations at low density are the conditions for which the linear 
theory is most valid but the Newtonian approach is least valid. Nevertheless, 
the simplicity of the Newtonian approach serves as a valuable guide in the sub- 
sequent treatment, and so far as the author is aware, it has not been previously 
developed in comoving coordinates. 


11 



For a collisionless fluid, such as a supergas of stars or galaxies, the formal 
approach is by way of the Vlasov equations , 30 as in plasma physics. This has 
been used by Gilbert , 29 and by Sweet 31 for counter-streaming fluids. Particles 
travelling an appreciable fraction of a wavelength in an oscillation period cause 
Landau damping. Since we are concerned with an initially structureless fluid we 
shall use the fluid approximation; this will tend to overestimate the rate of growth 
of perturbations. In this discussion the velocity components are everywhere single- 
valued, and only gravitational interactions are considered. 

Let the disturbed velocity, density, pressure, and gravitational potential be 

u-u + v, p -> p + S p, 

p p + Sp, p* p* + <p, 

where the small quantities are functions of r and t . In the usual coordinates t , 
rS/S 0 the linearized equations of motion and continuity are 

(— + u • V ,N ) v + v-V' u + V' ip + — V' Sp = 0, (15) 

\B t / P 

(JL + u-V' + V' -u^Sp + pV'-v = 0, (16) 

and terms quadratic in small quantities are neglected. The perturbation of 
Poission's equation is 
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The transformation to comoving coordinates is 


S o „ S o /. B . B . B \ 

— V = — [ l + + 1 , 1, 

S S \ 3r r d 0 rsin 6~d<pj 


— + u • V' = 
Bt 


3 

— + r 
Bt 


S_3 

3 Br 


d 

dt 


where the time derivative now follows the unperturbed motion of the fluid. In 
ordinary coordinates B/Bt commutes with V', whereas in comoving coordinates: 


_d 

dt 



With the use of comoving coordinates, (15) - (17) become 


— (— v] + V0 + 1 VSp = 0, (18) 

dt \S Q j p 

— + — v • V = 0, (19) 

dt \p ) S 

S 2 V 2 0 - 477GS 2 s p - 0, (20) 

where V • u = 3S /S Q , v • Vu = v S/S 0 , and pS 3 = constant. 



Taking the curl of (18), we have 


— (S a J)=0, (21) 

dt 


where V'/sv = { is the vorticity. Hence S 2 J is conserved. In the absence of 
pressure and density gradients both v and { are zero. The vorticity is therefore 
permanently zero and the motion is irrotational: 

v = (S 0 /S)V(p. (22) 

For the four unknowns 8p, 8 p,4> t and cp, there are the three equations (18) - (20). 
An equation of state is therefore necessary. Let pa p y , where y is the ratio of 
specific heats, then 


8p = y(p/p) Sp = c 2 Sp’ (23) 

and c s is the speed of sound. Taking the divergence of (18) and using (19), we 
find 




V 2 >// = 0. 


(24) 


It is now possible to write ^ = i/»(t) n(r,i9,0), and 

V 2 n+k 2 n = 0, (25) 
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where k 2 is the separation constant. With (5) and (7), Equation (24) becomes 

0+4 — 0 + ( — c 2 k 2 - 2 — +A^ 0 = 0, (26) 

s Vs 2 S 2 / 

and dots denote time derivatives. The remaining two equations are 

I il (S0) = 4 ttGp(P, (27) 

S dt 

S 2 k 2 0 = - 4 -ttGS 2 3 P- (28) 

The last three equations determine the fluctuations of density, velocity, and 
gravitational potential for an inviscid, irrotational fluid. They are identical with 
those obtained in Section 5.5 from the theory of general relativity for c 2 « c 2 . 

3.3 Normal Modes of Vibration 

An arbitrary disturbance in a scalar quantity consists of a superposition 
of normal modes given by (25). These modes can be constructed from plane 
waves 

n k oc exp i k • r 

The modes in spherical coordinates are given for comparison with those in Sec- 
tion 5.7 for curved space. 
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Let H - W(r)0((9)cl)(0); by separating the variables, Equation (25) in flat 


space becomes 


d 2 <t> 

d^2 


+ m 2 $ = 0, 


(29) 


1 d / . n d©\ ( , x m 2 

— sin 6 | + n(n + 1) 

sin# dfe 1 \ dd} V sin 2 # 


0-0, (30) 


J_ JL ( r 2 ^-\ + (k 2 - n(n + ^ V = 0. 

r 2 dr \ dr 


(31) 


The tesseral harmonics of n th degree and nth order are 


0 <t= Y^(#, <P) - (a nm e im * + b nm e- im ^) P" (cos #), 


for integral values of n and m , and the spherical surface harmonics of degree n 


are 


Y „(#,<£) = 2^ Y n <*>)• 

m-O 


(32) 


The solutions for the radial function 'i'(r), for n integral, are the spherical 
Bessel functions 
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(kr), «P n 



(kr). 


There are no boundary or periodic conditions to satisfy, and the only condition 
is that 'P must be finite everywhere. For kr 0, 


V n - (kr) n ./ [1.3.5. ...(2n + 1)], 


- [1.3.5. . .(2n - 1)] /(kr) n + 1 . 


and kr -oo, 


*P n - (kr)' 1 cos |kr - i (n + 1) v , 

- ( kr )’ 1 si " J kr ( r * + ') " j * 


and hence V is rejected since it diverges as kr -0. The radial functions form 
a continuous set having all eigenvalues of k > 0 for each value of n . In particular, 
for n = 1, 2, and 3: 


*P Q = (kr)' 1 & 

= (kr)' 1 [(kr)' 1 & _C] , ( 33 ) 

V 2 = (kr)' 1 [ (3(kr)' 2 - 1} & - 3(kr) _1 C], 
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sin kr, C = cos kr. The spatial disturbances are represented by summations 
and integrations over the complete set of wave functions 

n = \ (kr) Y” ( 0 , 4 >). ( 34 ) 

4. GRAVITATIONAL INSTABILITY 
4.1 General Criteria 
Let 


= S m Spa Sp/p m ^ 3 . (35) 

The existence of time growing modes of an arbitrary quantity such as a m does 
not necessarily imply instability in a nonstatic universe. For example, if a = 8p 
grows, an expanding universe is unstable, but a contracting universe is stable 
if p increases more rapidly than Sp. The only time growing quantity that denotes 
unambiguously a changing configuration is the contrast density a 3 = S p/p . Sup- 
pose that a m grows in time; then in an expanding (contracting) universe m = 3 is 
necessary and m ^ 3 (m 3) is sufficient for instability. Thus the growth of the 
gravitational potential (27) = a a 2 , is sufficient only for instability in an ex- 
panding universe. 

When the growth time of a mode is greater than the age of the universe it 
cannot contribute to a significant change in configuration. For instability we 
require 

a 3 /a 3 » |S/S|, (36) 
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oc 


or, if a ac S ±C (+1 for S > 0, - lior S < 0), 


f » 1. 


Even then a clear case of instability requires an adequate amplitude of the initial 
disturbance. For arbitrarily small initial conditions a clear case of incipient 
instability requires a 3 a exp a, t, a> » js /S | • 

The Newtonian equations of the unperturbed cosmological models, and the 
linearized Newtonian equations of the perturbed models, are identical with the 
corresponding general relativity equations for small pressure. The Newtonian 
equations are limited to flat space and cannot be used to determine the behaviour 
of modes in curved space. This is not a severe limitation as the eigenvalues k 
in (26) are multiplied by the speed of sound c s , and the time dependence of a m 
can be discussed for a given value of kc s , 

From (26) the general equation for a m is 


a +(8-2m)_a + 

m v ' _ m 


K 2 + (6 - m) £ + (4 - m) (3 - m) — - A 
S S 2 


a m “ 0 , (38) 


where K = c s k S 0 /S. Equation (26) is for m = 2; and for m = 3: 


oi 3 + 2 - a 3 + (K 2 - 4wGp) « 3 = 0. 

S 


(39) 
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This equation has been derived by Bonnor 24 and van Albada 25 for radial type 
perturbations and has been discussed more generally by Savedoff and Vila. 27 
From (27) the velocity potential is given by 

<x 3 + cpk 2 S 2 /S 2 = 0. (40) 

For m = 4, 

a ’ 4 + ^K 2 - — 77 G p - — a 4 = 0. (41) 

Unless stated otherwise we shall assume that the cosmological term A is 
zero. When 

K 2 > 8ttG/V3, ( 42 ) 

a 4 is oscillatory, and hence a m is also oscillatory for all m. For k = 0 (zero 
energy) this condition is e 2 > (Sft/ S 0 ) 2 , where ft. = k' 1 , and the velocity of sound 
exceeds the expansion or contraction velocity (1) at distance ft . 

Since pa p?, it follows 

'I =<o (S 0 /S) 3r “ 3 < 43 > 

or K 2 = K 2 (S 0 /S) 3T 1 where K Q = c sQ k at S = S Q . Equation (41) is now 

a 4 + [K 2 (S 0 /S) 3r ' 4 - /3/S 3 0 ] (S Q /S) 3 a 4 = 0. (44) 
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This eauation is useful for discussing qualitatively the behaviour of a m as a 
function of time. So far S 0 is arbitrary; let S = K 2 at S 0 = P/Sff °r 

c^ 2 = 8ttGp 0 /3, (45) 


and (44) becomes 


a; + [(S 0 /S) 3r_4 -l](/3/S) 3 a 4 = 0. (46) 

Hence, a 4 has marginal stability at S = S 0 - The results are shown in Table 1. 

In an expanding (contracting) universe a m is a growing function of time for m > 4 
(m < 4) for all y. More precise conditions are given in Section 4.3 from the 
solutions of (38). At this stage we are unable to deduce any conditions concern- 
ing the growth of the contrast density (m = 3) in the important case of an expand- 

inrr nmvorco 

*■^^5 U.4JUL* ^4. W • 

A marginal state for any value of m is obtained by using 


e 


1 


(S/S 0 ) 2m - 8 dt, 


in (38): 


^ am + K 2 + (6 - m) + (m - 4)(m - 3) ■?_ 

de 2 L S S 2 _ 


. 2 

€ a - 0 • 

m 


For m = 3, this is van Albada's 25 equation 


d 2 a, / S oV 

+ (K 2 - 4 ttG p) [ — — ) a =0, 

de 2 VS / 


(47) 
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TABLE 1 


Time Dependence of a m 


m < 4 

C S > 0 ~ l, (?) 

L s < 0 ~ T) (T) 


m = 4 m > 4 

(T) ~t, (t) 

~, ( t) ~ i. (?) 


7 = 


4 

3 


4 

7 > — 

3 



?, (?) 

t , (T) 

?. (~ t ) 

t , (~ T ) 


t, (T) 
t, (t) 

T . (~) 
T. (-) 


T > (t) 

?. (?) 

T, (~ T) 

?. (- i) 


The symbols denote: S > 0, expansion; S < 0, contraction; t , growth; i , 
decay; oscillating; without brackets, S > S 0 ; with brackets, s <S 0 * For 
example: (~i) means decaying oscillation for S < S 0 . 
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and therefore, Jeans' criterion of marginal 'stability': K 2 = 4ttG p, holds for 
a 3 ( £ ) • However, e is not a linear function of time, and in general this marginal 
state does not give an unequivocal stability criterion. 

4.2 Jeans' Criterion 

Jeans’ treatment 21 of gravitational instability resembles Lord Rayleigh's 
formulation of the problem of oscillations in a fluid of positive and negative 
charges. 32 Jeans assumes that a uniform, gravitating fluid is in a stationary 
state (as for a neutral plasma) and S = 0. Hence (39) becomes 

a" 3 + (c 2 k 2 - 4ttG p) a- 3 - 0, (48) 

and S = S Q . For a 3 a exp icot this equation gives the dispersion relation 

a? - c 2 k 2 + a) 2 . (49) 

A similar relation holds for electrostatic oscillations in a plasma, and w =(47me 2 /m) 1/2 
is the plasma frequency for electrons of number density n and charge to mass 
ratio e/m. In Jeans' dispersion relation the 'gravitational frequency' is imaginary: 
oj p = (_4 7 tG p) 1/2 * There is thus a marginal state 


k - (4 ttGp ,/c 2 ) 1/2 > (50) 

J 5 

and for k > kj , q> is real and the disturbance oscillates at constant amplitude; 
and for k < kj , co is imaginary and the disturbance grows exponentially in time. 
Jeans' stability criterion is therefore k > kj, or ^ where X. = k* 1 , ^ = k^ 1 . 
The Debye length Xp = plays a similar role in plasma physics and 
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disturbances of k <\ in many cases tend to be stable. In a Jeans’ sphere (as 

compared with a Debye sphere) of radius k^ the sum of the thermal and potential 

energies is of the order e 2 pk* - 4nG p 2 k*=0. In a sphere of radius k < k , the 

S J J J 

thermal energy predominates and collective interactions are of little consequence. 

Jeans' analysis suffers from the defect that in general there is no initial 

stationary state in a uniform nonrotating fluid. When k «k J9 then a? » S 2 /S 2 , 

and the dispersion relation (49) is an acceptable approximation. But as k increases, 

the oscillation period also increases and is infinite at k - k ; and when k > k , the 

J ^ 

e-folding time i or 1 of a disturbance is comparable with the collapse time i(2/3) 1/2 af 1 
of the system. Thus, in the range of interest: k > k } , the dispersion relation 
(49) fails and we must fall back on solving (39). It is seen that Jeans' instability 
criterion: k > kj , is necessary for expansion (S > 0) and is sufficient for con- 
traction (S < 0), but in neither case is it both necessary and sufficient. 

In Einstein's static universe (which has its analogy in Newtonian theory, 2 
provided p « pc 2 ): S = 0, S = 0, at S = S Q , and according to (5) and (7) this is 
possible for k - 1 and 


A = S 0 2 = 4 7T G p 0 . 


(51) 


Equation (38) for a m becomes 

A" + (c 2 k 2 - A )p. = 0, (®2) 

where p. stands for a , i/< , or <p. As Bonnor 24 has shown, Jeans' dispersion re- 
lation (49) holds true without modification in Einstein's static universe. The 
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I 


cosmological term neutralizes the gravitational field and there is a similarity 
with the neutral plasma state. If the pressure is zero, as in a fluid consisting of 
dust particles having no peculiar motion of their own, the velocity of sound is 
zero and all modes grow as exp A 1/2 t. 

4.3 General Equations 

To solve (38) we use the relations (11) and (12) and take X as the independent 
variable. For k = 0 (zero energy). 


d2a m + (14 - 4m) 
dy 2 x dy 


+ 4 


L(y) + 


(m - 2)(m - 9/2) 


a - 0 , 


L(y) - 




6 


(53) 


k - 1 (negative energy): 


d 2 a da 

EL + (14 _ 4m)cot y — 

dy 2 dy 


L(sin y) - (m - 4)(m - 3) + 


(m - 2)(m - 9/2) 


0 , 


sin' 


L(sin y) = K 2 /3 2 (sin y ) 
u o 


6r 


2 / . 67 
(sin y) 


(54) 


and for K = -1 (positive energy): 


d 2 < 


da 


+ (14 - 4m) co thy + 4 

y dy 


d - 2 


L(sinhy) +(m -4)(m - 3) + 


(m - 2) (m - 9/2) 


S 1 


inh 2 


a - 0, 


L(sinhy) - K 2 f? (sinh y^ ) 6r '^/(s inh y) 


67-6 


(55) 
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These equations are used to consider the growth rates in (i) a cold universe, 
(ii) an isothermal universe (y = 1), and in the cases when (iii) 7= 4/3, and (iv) 

7 = 5/3. 

4.4 Cold Universe 

If the pressure is zero, as in a universe containing dust particles having no 
peculiar motion of their own, all time growing perturbations have their maximum 
possible rate of growth. A cold universe is the most unstable of all, and is not 
only a very simple but is also a very interesting model to study. In this model 
K 0 = 0 in (53)-(55). 

For k = 0 (zero energy) the solution of (53) is 

a m - A S m ' 9/2 + B S'"' 2 , ( 56 ) 

A, B are constants. Hence, during expansion (S > 0),a 2 is either constant or 
decays, and during contraction (S < 0),a g/2 is either constant or decays. For 
S > 0 (S < 0), a grows for m > 2 (m < 9/2). Using only the growing terms 

m 

S > 0: 8p/p oc S a p " 1 ' 3 , 

S < 0: Sp/pocS‘ 3/2 a p 1/2 , 

and these growths are the same for all modes. 

The absence of exponential growth is typical of all non-static models. Ac- 

3 

cording to (36) neither -l = 1 ( expansion) nor -fc = — (contraction) satisfy the 
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condition (37): -t » 1, for a convincing case of instability. Because the universe 
is expanding the result *t = 1 is of particular interest. It shows that during the 
life-span of the universe in which the Newtonian theory is valid small initial 
disturbances tend to remain small. For example, if initially p l ~ 1 g cm' 3 and 
^P l / p x ~ 10 ‘ 15 , then at p ~ 10 30 gem' 3 it follows that Sp/, 0 ^ 10 ' 5 . 

In the k = 1 (negative energy) or oscillating model the solution of (54) is 


a m = ( sin *) 2m ' 7 [AP» (i cot y) + B Q‘ (i cot y) ] , 


(57) 


where 0 £ X 1 tt, A and B are different constants, and p£, are the associated 
Lengendre functions. Since P 1 (-ix)= -P 1 (ix), Q 1 (-ix)= Q 1 (ix), we need 
consider only the range 0 £ y £ -j tt.. Using the relations 33 


Pj (ix) = - 3(1 + x 2 ) 1/2 x, 


Q 1 (ix) = i (1 + x 2 ) 1/2 [3x cot' 1 x - (3x 2 + 2)/(l + x 2 )] , 


(57) becomes 


a m = & 2m ' 9 [AC + B(3A - & 3 - 3 y C)] , 


(58) 


& = sin y , C = cos y. As y - 0, this equation is identical with (56). 
In the k = -1 (positive energy) model the solution of (55) is 


a m = (sinh y) 2m ' 7 [AP 1 (coth y) + BQ 1 (coth y)]. (59) 
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From the expressions (x > 1) 


$ 


(x) = 3(x 2 - 1) 1/2 x, 


(x) = (x 2 - 1) 1/2 


3 

2 


x ■in 


x + 1 
x - 1 


3x 2 - 2 
x 2 - 1 


it is found 


a m = S 2m - 9 [AC + B(3S + S 3 - 3yC)L (60) 

S = sinh y , C = cosh v As y - 0, (56) is recovered; however, as y -oo, a m = 

A s 2m ' 8 + b s 2m ' 6 > and in an expanding (contracting) universe a m grows for m > 3 
(m < 4). 

The solutions (56), (58), and (60), for a 3 = 8 p/p, are shown in Figure 1 with 
the constants A and B equal to unity. In an expanding universe the growth is 
least, as one would expect, in the case of positive energy of k = - 1 . The rate 
of growth is greatest for k = 1 , but the growth is limited because the universe 
oscillates. 

4.5 Isothermal Universe 

Under certain conditions it is possible for a fluid to expand and contract 

ryry 

isothermally. Furthermore, it may be possible for regions of increased 
density to be cooled by radiation, thus tending to make the fluctuations isothermal. 
When y = 1, the L's in (53) -(55) equal (c s kS 0 ) 2 . At S = S Q , let 

c 2 0 k 2 = 8 ttGp 0 /3 (61) 
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AMPLITUDE 





Then, for * = 0, (c s0 k S 0 ) 2 = fi/S 0 , and 


= S 


m* 1 3/4 


[ A J 1 (2S 1/2 /s 1 0 /2 ) + B J l(2 S l/2 /s\ /2 )j 


(62) 


and therefore 


S « s 0 : a m = A S"' 972 + B S"" 2 

S » s 0 : a m = S"" 772 [A sin 2(S/S 0 ) 172 + B cos 2(S/S 0 ) 172 ]. 

At high density or long wavelengths a behaves as in the cold universe, and at 

m 

low density or short wavelengths a ? , 2 oscillates. 

For k = 1, 

a m = ( sin y) 2ro ' 1372 [AP^ /2 (cos Y ) + BQ 5 / 2 (cos y) ] (63) 

and k = -1, 


a m = (sinh y )2-"-i3/2 [AP^ /2 (cosh y ) + BQy 72 (cosh y)] (64) 

where X. = - 1/2 ±(1 + 4 k c 2 q k 2 S 2 ) 172 = - 1/2 ± (5 + 4 k S 2 ) t/2 and S 0 occurs at 
S Q . When 2c s0 kS 0 < 1, these equations are similar to those of the cold model, 
and when 2c s0 kS 0 > 1 they resemble the zero curvature solution (62). Savedoff 
and Vila 27 have derived results similar to (63) and (64) in terms of hypergeo- 
metric functions. 
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4.6 Arbitrary y 

The general solution of (53) for K - 0 and y arbitrary is 


a m a S"' 1374 J ±p (2 [ S / S 0 ] 2 ' 3y/2 ), 


( 65 ) 


P = 5/2(4 - 3 y), 


where (61) holds at S = S Q ; and for y = 4/3, 


a a S m± ^' 13/4 , 


( 66 ) 


, 3c; o k^ 

\16 8ttG p 0 


1/2 


and c s is the speed of sound at density p 0 . From (65) 


4_ 

T 




y < — i S « S Q : 


> a = AS " 1 ' 972 + BS"’ 2 , 

/ m 


y > — . S » s 0 : 

3 J 


(67) 


y < , S » S Q : a m = s m ' n ' 1374 (A sin x + B cos x), 


y > — , S « S Q : a m = s nl+n ‘ 13/4 (A sin x + B cos x), 

3 


( 68 ) 
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where n = 1 - 3/4 y, x ~ 2(S/S 0 ) 2 ' 3y/ 2 . The general solution has been given by 
Savedoff and Vila. 27 

A photon, neutrino, or relativistic gas has a ratio of specific heats of 4/3. 
When the pressure of such gases is dominant but their density is small com- 
pared with the total density of the fluid, the fluid as a whole has y = 4/3 and can 
be treated with the Newtonian approximation. When ( k = 0) 

c so k 2 > — 7rGp 0 , (69) 

6 

a, oscillates: 

13/4 


a m « S’ 1 ' 13/4 e ±l s , 


(70) 


where i /jl = X in (66). In the case of k = ±1 


a m ~ 7/2 [ AP^(i cot K y) + BQ^ (i cot K y) ] , (71) 


and 


\ _ i. + (25/4 - 3c 2 q k 2 /277G/o 0 ) 1/2 . 

, Q^are nonperiodic for condition (69). 

When y = 5/3, /< = 0, (65) is 


a a 
m 


13/4 


J± 5/2 ( 2S Y 2 /S 1/2 > 


(72) 
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as shown by van Albada. 2 5 There appears to be no simple analytic solutions for 
k - ±1 for arbitrary y , nor is there in the important case of y = 5/3. 

5. RELATIVISTIC COSMOLOGY 
5.1 Unperturbed Models 

In its unperturbed state we assume that the universe is homogeneous and 
isotropic, and the metric is given by the Robertson-Walker line element 

ds 2 = dt 2 5! (dr 2 + r 2 d Q 2 ), (73) 

c ! ( 1 + 5" <r2 ) 2 

where dfi 2 = d8 2 + sin 2 0 d cj? , r, 9, 4> are comoving coordinates, R(t) has the 
dimensions of length, and k = 0, ±1 is the curvature constant. The energy- 
momentum tensor of a perfect fluid is 

Tj = (pc 2 + p) g k . u k u 1 — SI p , (74) 

in which pc 2 is the energy density, u 1 the four velocity, and p is the isotropic 
pressure. For a fluid that is stationary in the comoving system, u° = 1, 
u 1 = u 2 = u 3 = 0, and the components of the energy-momentum tensor are 


Tq = ^ c2 ’ T 1 = T 2 = T 3 = - P- 


(75) 
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All that now remains is to determine R(t) in (73) with the Einstein equation 

R- - — 8! + S* A = - ( 877 G/C 2 ) T 1 . (76) 

j 2 * ^ J j 

In this equation, Rj is the contracted Riemann- Christoff el or Ricci tensor (its 

t 

further contraction R^ is shown explicitly to avoid confusion with R(t) ), and the 
cosmological term A is included as in the Newtonian treatment for historical 
interest. 

Equation (76) is readily solved using the line-element (73) and the com- 
ponents (75) of the energy-momentum tensor. 1,34 The following method however 
is particularly simple. Transforming to the coordinates: 

x = r sin 6 cos 0, 

y - r sin 8 sin 4 >, 

z = r cos < 9 , 

r 2 = x 2 + y 2 + z 2 , the line element (73) in the immediate neighborhood of r = 0 


becomes 



At any desired point r = 0 the coordinates are natural, 35 as in the Minkowski 
metric. For these coordinates the Christoff el symbols 


i / p S B e qs 3gp C 

r: =-g rs f— + 


pq 


BxP 9x s 


are 5 * 


r s. = R * 2 a = k <»/». 


r“ = - rt = -r a =Iix‘ 

00 0a aa 


(78) 


As r - 0, the only surviving Christoffel symbols are T“ a and r£ a , and the spatial 
derivatives of (78) are 


B 

Bx a 




_a_ r * = JL 

3x a aa 2 


(79) 


Since 


= 


Bx' 


r 1 , - -iL- r,^. + r m o r^„ - rr 


k'l 


Bx 1 


kj 


k^ V 


kj 




(80) 


*The convention adopted is that Latin indices assume all values 0,1, 2,3: Greek indices \, fi, v 
assume only the values 1,2,3, and a, (i , y are used when there is no summation and a y fi yy . 
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we find from r| = g 11 ^. that 


Rjj = 3 R/R, 

R“ = (RR + 2 R 2 + 2 k c 2 )/R 2 , (81) 


and all other components vanish. From (76) we now obtain the well-known 
equations 


R 2 = i ( 8 ttG /0 +A) R 2 - k c 2 , 
3 


(82) 


i. (R 3 P) + L =0. (83) 

dt c 2 dt 

Every observer can adjust the origin of his comoving coordinate system to give 
r = 0, and therefore (82) - (83) apply to all comoving observers. Equations (82)- 
(83) are identical with (6)-(7), provided R = cS and p/c 2 is vanishingly small. In 
the relativistic treatment we shall use R rather than S. 

5.2 Equation of State 

We have two equations, (82) and (83), for R , p, and p, and therefore require 
an equation of state: 


p = (v - 1 ) pc 2 . 


(84) 
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In any range of density in which v is constant, (83) becomes 


If, now 


(82) can be written as 


For the three models k = 


K = 


K - 


K = 


where 


/OR 3v = constant. 


P v = STTGpR^/a, 


(85) 


R 2 



A R 2 - k c 2 . 


( 86 ) 


0, 1, -1, and A = 0, it is found that 


0: R = (fi v c' 2 y 2 ) l/(31 ' . 

1: R = (j3 v c' 2 sin 2 y) 1 ^ 3 v ' 2 \ (87) 


- 1 : R = (fi v c 2 s in h 2 y) 1 ^ 3 "' 2 \ 


dt 2 R 
3 v - 2 F 


( 88 ) 
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In the relativistic treatment the pressure contributes to the mass density of 
the fluid. If the kinetic and non-gravitational interaction energies of the con- 
stituent particles of the fluid are small compared with the particle rest masses, 
as in the universe in its present state, there is justification for assuming v - 1 
as in the Newtonian theory. However, there is no justification for assuming that 
the pressure gradients in a perturbed universe of p « pc 2 are also of negligible 
effect. If, at low density, instead of (84) we adopt 


Sp = (v - 1) c 2 Sp 


(89) 


where Sp = c 2 8 p, then 


^=l+c 2 /c 2 = l+ yp/pc 2 , (90) 

and c 2 « c 2 . In an isotropic photon or neutrino gas, or in an isotropic fluid in 

which the particles and their fields have energies large compared with the rest 

masses, v attains its maximum value 36 of 4/3 and p = — pc 2 . The range of v 

3 

is therefore 1 £ v < 4/3. (Zel'dovitch 37 has proposed increasing the upper limit 
to 2, but it is possible that this is unrealistic. 38 McCrea 39 suggests a lower 
limit of -2 as a physical explanation of a steady state universe, but this is out- 
side the scope of our discussion.) 

In the case of a gas containing degenerate electrons, the electrons are rela- 
tivistic when their Fermi energy E e exceeds m e c 2 (m e is the electron mass). But 40 


38 



( 91 ) 



where a e is the mean interelectron distance: — 77a 3 n e = 1, n e the electron density, 
and K e - h/m e c the Compton wavelength. Hence, E e > m e c 2 occurs in hydrogen 
at p > 10 6 g cm' 3 , approximately. It is reasonable to suppose as a first approxi- 
mation 41 that the Fermi energy level is comparable with kT, where T is the 
temperature (in this section k is the Boltzmann constant). As the density rises, 
electromagnetic and weak interaction populate the photon and neutrino states, 42 
thus ensuring that the Fermi and thermal energies (after allowance for rest 
masses) of the fermions and bosons are of comparable magnitude. The leptons 
and photons contribute most of the pressure and yet contribute only a small 
fraction of the total energy pc 2 of the fluid. Therefore (90) is used with y - 4/3. 
The lepton and photon energies become significant as the Fermi energy approaches 
m n c 2 (m n is the nucleon mass), and according to (91) this occurs for a e ^ \ , where 
k n =fi/m n c (ra n is the nucleon mass), or when p -v 10 16 g cm" 3 . (The condition 
E e ~ kT limits electron pair production and the number of electrons does not 
greatly differ from the number of nucleons. Inverse /3 decay: 40 - 43 p + e" - n + v , 
which is of importance in stellar physics, plays only a minor role because of the 
dense neutrino-antineutrino background. The results of other studies, 44 in which 
it is assumed that the temperature is low, are also inapplicable.) The density of 
nuclear matter is p n ^ 10 15 g cm* 3 . For p < p n the pressure in general is small 
compared with the energy density and v has a value close to unity. When P ^ p n 
the composition of the fluid is complex, and at higher densities progressively 
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increases in complexity. However, whatever the nature of the fluid at super- 
nuclear density, it seems likely 38 that v will have a value close to 4/3. By these 
crude arguments it is possible to distinguish two states of the universe: 

subnuclear state: p < p , v ~ 1, 

. . . 4 

supernuclear state: p > p n , v ^ , 

and in each of these states v is roughly constant. The principal change in v 
occurs in the transition through nuclear density. 

Neglecting the transition region of v there are now altogether six models: 

Subnuclear-models . v = 1. These have been previously given in (11)-(12), 
but are now repeated using (87)-(88): 

k - 0: R = /3j c' 2 y 2 , 

k - 1: R = c -2 sin 2 y, 

k - - 1: R = c' 2 sin h 2 y. 


and dt = 2Rdy/c. 

Supernuclear models, v - 4/3. Again, from (87)-(88): 
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* = 0: R " ^4/3 c' 1 y. 

* = 1: R = y8j/J c-1 sin y, (93) 

K = - 1: R = ^ 4/3 c ' 1 sinh y, 


and dt = Rdy/c. 


The value of y at the density p is estimated in the following way. At 

n n 

present in the universe ~ R Q c 2 , and hence B ~ R n ~ R n R o ° 2 • Also, 
y ~ R n c /p\% and therefore 

y ^ (R n /R 0 ) 1/2 ~ (P 0 /P n ) 1/6 . (94) 

n 

From the present density of p Q ^ 10' 30 g cm' 3 it follows that y ^ 10' 8 . The 

n 

supernuclear models (93), in which y < y , have therefore negligible difference 

n 

and for simplicity we can assume k - 0. 

5.3 Linearized Equations of Perturbed Models 

We consider small departures from the metric (73) as the result of displace- 
ments of the fluid. A perturbation treatment of the cosmological models, as 
distinct from a static and flat metric, 34,45 encounters the slight complication of 
non-vanishing Christoff el symbols. The most general treatment has been given 
by Lifshitz 46 and includes rotational motions and gravitational waves. Owing to 
the conditions imposed on the metric tensor, simple ir rotational motions were 
excluded. The Newtonian approach shows us that in the simplest treatment 
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these are the motions of most inerest. A treatment of the problem, analogous 
to the Newtonian treatment of a uniform fluid, does not appear to have been pre- 
viously attempted in a simple and straightfoward manner. 

Small variations in the metric tensor are expressed as g + Sg. k , where 
the g jk are given by (73) and 

Sg jk =h jk - (95) 

We assume h. k and its derivatives are everywhere small, and that quadratic and 
higher order terms in small quantities are negligible. Thus the unperturbed 
tensors g jk , g ik are used for lowering and raising the indices of h lk ,h jk : 
hi = g ik h. k = g. k h ik , and in effect h jk is a tensor field in the unperturbed g jk 
space. Since g.. g ik = Si , to a first order 

J K J 


S(g jk g lk ) = hi + g jk Sg ik = 0, 

and therefore 

Sg ik = _ h ik . (96) 

It is more convenient to use Einstein's equation (76) in the alternative form 

Rj = - (8 77G/C 2 ) (t; - \ S.t) + s. A, (97) 


and its perturbation is 


SRj = - (8rrG/c 2 )S ^Tj -- 8} T) 


(98) 
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The perturbation S R‘ is evaluated in terms of hj as follows. From (80), 


5R « = sr ^ ‘ r ^‘ sr <- " r l 8r ^ - r i 8r :-e) 

- {-ti 

thus giving the Palatini equation 47 

S*J k = Sr U;, - Sr Li’ <") 

where a semicolon denotes covariant differentiation. Also, to a first order, 


sr; = S(g ik R ik ) =6 ,k SR ik -i<R5. 


(ICO) 


and with (99), we obtain 


g ik c 8 r’u-fc s j - 


A 


- SR j +K#r 


(101) 


The perturbed Christoffel symbols are 



Hg" 1 ' g rl r 


kj 


) = g^ r S(g r 


r s ) -r* 

s kj ' kj 



s 


= i /, (fv + 

2 8 \3x j 3x k 



- n- 

kj 



43 



thus giving Lifshitz’s 46 equation 


"t 4 S' ♦'WW' 


( 102 ) 


't 1 't 

In particular, 8 = -h k , where h is the trace h^. 


From (101)-(102) we have a linearized differential equation in hj: 


6‘ k (N„, - •^u- h W*«** h U.= l < s,, > *><*?>• 


4 




(103) 


The g.. and R? are known, and 8R1 is given by (98) in terms of the perturbed 
energy-momentum tensor. 

The trace of t| is T = pc 2 - 3p, and therefore 

8T = c 2 Sp - 3 Sp. (104) 


Furthermore, g k . u k u j = 1, and from §(g kj u k u j ) = °» 

h kj u k u j + g kj 8u k u j + g kj u k Su J = 0 

where u° = 1, u“ = 0, and therefore 



Also, g^i^iF = S' 8®, and 


8 (gkj ^u 1 ) 


= Si o h ° 


S°u„ 


)U + 


gkj 8 u k 8^ u c 


We therefore find, from (74) 


ST' = S (pc 2 + p) gk . u k u' + (pc 2 + p) 8 (g k j u k u' ) - S' Sp 


= si 


s° S (pc 2 + p) 


- S' 


Sp + (pc 2 +p) (S' h® 


S° 

J 


u 0 Su« 


gkj 8 <J k 8 ju°). (106) 


From (105) and (106) it follows that the components of ST' are 


ST« = c 2 Sp, 


ST® = (pc 2 + p) (h® + g aa Su a ) , 


(107) 


ST“ 

a 


-Sp, 


STJ = 0, 


where a,/3 = 1 } 2, or 3, a / yS, and no summation. From (98), (104), and (107), we 
now have 
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SRO = _ 47tG S(p + 3p/c 2 ), 


SR° = -8ttG(/ 3 + p/c 2 ) (h° + g aa Su a ), (108) 

8R“ = 477GS(p - p/c 2 ), 

and 8R^ = 0. In addition, the components of R| are 

RO = -4ttG Go + 3 p/c 2 ) + A, 

R“ = 4ttG (p - p/c 2 ) + A. (109) 

Collecting together the equations (103), (108), and (109), we have, for i = j = 0: 
h ; oo - 2h \-,ol + h 0 ;^m = " 877G ( h o + S ) (P + 3p/c 2 ) + 2h° A; (110) 

for i = 0, j = a, (a. - 1, 2, or 3; no summation over a): 

h - s - s + h° 9 ~ 

;0<x u 0;a/ ,1 a;0 / t T ® a ; Am 

- 8?7G [h° (p + 3p/c 2 ) + 2 (p + p/c 2 ) g aa Su a ] + 2 h° A; (HI) 
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for i = a, j = ol: 

g aa (h ;aa -2hf. a ^) + g^ m h“.^ m = 877G(h“ + S) (p -p/c 2 ) + 2 h“ A; (112) 


and for i = a , j = /3 , (a^/S): 



Altogether, (110)-(113) provide ten equations for the determination of the ten 
unknowns: h k . (four component's can be discarded by coordinate transformations), 
Su a (since Su° = - 1/2 h°), and Sp (Sp is given by an equation of state). 

5.4 Irrotational Motion 

For irrotational motion it can be shown that (110)-(113) reduce to three 
equations determining </», cp, Sp, as in the Newtonian treatment. 

The simplest procedure is to transform to local coordinates where the 
Christoffel symbols are given by (78). In addition, we adopt a system of coordi- 
nates in which 


h 


Oa 


= 0 , 


(114) 


and 



(115) 
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The only equation involving the velocity 8 u a is (111). This equation now becomes 


3 

3x“ 




30 


3t 


- = - 16 G (p + p/c 2 ) v a , (116) 


where g aa Su a = v a . For irrotational motion v a;/3 = vg ;a , or 


v = 3cp/c 2 3x a , 


(117) 


and this is possible when the left-hand side of (116) is a derivative with respect 
to x a . Hence, 0 a = 30/3x a ; and since 0can be absorbed into cp, we assume 0 = 0. 

Equation (113) is now 


^|2l?R 2 +R 3 ^| +Rc 2 (3/<-V 2 ) 


= Aj 




For (111) we take i = 0, j = a, and subtract i = 0, j = fiz 


8-O0.-p/c>) t 2A + | + E t f^-n*_ 


(hj - h 'b -- 


ts\ (h; - h|) = o, 
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where □ 2 is the d’Alembertian operator. From A 2 = 0, A 2 (h“ - h^) = 0, it is 
seen that h”. and h“ - h| are propogated independent of any fluid perturbation. 
Since hg and h£ - h| are zero in the absence of perturbations, and we require 
only those hj which result from fluid disturbances, we can choose coordinates 
in which = 0, h“ = h^. Thus, all nondiagonal components of hj are zero, and 

h° = - h} = - h 2 = = 2 0/c 2 . (118) 

The factor 2 c' 2 is chosen so that \f> is equivalent to the gravitational potential of 
Newtonian theory. For small irrotational fluid disturbances the line element is 
therefore 


ds' 




dt 2 - 


R 2 (1 - 20c' 2 ) 

9 / . 1 9 \ ^ 

c* i + — k r* i 

V 4 ) 


(dr 2 + r 2 d.Q 2 ) 


(119) 


A transformation of the coordinates x* -> x' * = x 1 + e 1 , where e* are small 
quantities, leads to g jk + yj k , where 


^jk ^jk e j. k € k; j * 
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The line element (119) is unchanged with the Killing equation 


e . +£...- 0. 

J ; k k; j 
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From 


r|vi clx^ 

(gjk + h ,k)“77 TTT = (8 jk + y jk 


. dx' i dx' k _ j 
} ds ds 


it is found de°/ds = de°/dt = 0. Also, for Su a to remain unchanged, de a /dt= 0. 
Clearly, these conditions are not limited to infinitesimal transformations, and 
any x 1 - x' 1 leading to 0 (x 1 ) - 0 ' (x' *) is admissable 

5.5 Equations of Irrotational Motion 

With the line element (119), the only surviving equations of (110)-(113) are 
40 + 12RR' 1 0 - O 2 0 = 4ttG (20 + c 2 8 ) (p + 3p/c 2 ) - 20A, (120) 

4RR' 1 0+4(RR‘ 1 + 2R 2 R' 2 )0 + □ 2 0=47tG(2 0-c 2 8)(/O- p/c 2 ) + 2 0 A, (121) 

— (R0) = 477 G(P + p/c 2 ) cp. (122) 

R dt 

These equations are obtained either by working through the covariant differentia- 
tions, or more simply using local coordinates. The results have been checked 
with Dingle's 49 formulae for an orthogonal line element. With 

V 2 0 = - k 2 0, 
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(120)-(121) become 


□ 2 ^ = 







4ttG c 2 Sp=-3 RR' 1 \p - (k 2 c 2 + 2 AR 2 + 3R 2 - 3/<c 2 ) R' 2 0, (123) 


4 77 G Sp = 0 + 4RR' 1 0 + (2RR + R 2 - *c 2 ) R' 2 0. (124) 


The fluid motions are therefore parallel to the pressure and density gradients. 
With p - (v~l) pc 2 , for constant v, we obtain from (123)-(124) 

0 + (1 + 3t ) RR' 1 0 + [(y- l)(k 2 c 2 + 2 AR 2 ) + 2 RR 

+ (3v-2)(R 2 -kc 2 )] R' 2 0=O. (125) 

At subnuclear density v - 1 + c 2 /c 2 , c 2 « c 2 , and therefore 

7 a k /’ 

0 + 4 — 0 + 

R r 


- 2 — + A 0 = 0, 

_\ R / R2 


(126) 


since 2 R R + R 2 = AR 2 - k c 2 . This equation is identical with the Newtonian 
result (26) with the transformations R -> cS and k 2 - k 2 S 2 /e 2 . Because Sp = c 2 Sp, 
it follows from (124) and (126) that 


51 



k 2 0 = - V 2 0 = - 477 GR 2 Sp, 


(127) 


which is Poisson’s equation. Thus, (122) with p = 0, (126) and (127) are the 
Newtonian equations (26)-(27). 

At supernuclear density v = 4/3, and from (125) 

0 + S^\fi + (^- - 4 + 2 a) 0=0. (128) 

R \3R 2 R 2 ) 

This equation is discussed in Section 5.8. 

5.6 Previous Work 

The problem of perturbed cosmological models has on the whole received 
scant attention in the theory of general relativity. The Newtonian approach to 
the problem is more popular in spite of its limitations. In the case of irrotational 
motion the two approaches lead however to identical results for relatively low 
pressure. 

Lanczos 50 derives an equation similar to (103) and adopts a generalized 
de Donder 51 coordinate condition 


(h 1 S 1 h) = 0, 

\ 1 2 > );i 


(128) 
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which is resembles (115) with 9 a = 0, jna. However, a scalar wave equation of 
zero nondiagonal components hj, and h° = -h“ = - -Lh , is impossible according 
to (128). The coordinate condition (128), commonly adopted in the study of gravi- 
tational waves, 34,45 yields in general tensor wave equations. 

Irvine 52 uses Lanczos* equation and shows that for irregularities of scale 
length ft , such that ft « R , ft R/R « c , an approximate form of (103) is 

□ 2 hj = 2 S Rj = -16 tt G 8 ^Tj - i- Sj T j . (129) 

Also, when the fluid velocity is small compared with the velocity of light the 
principle component of (129) is 


□ 2 h° = -SttGSCp + 3p/ c 2 X (130) 

In this approximation 8 p/p, Sp/p, h°/ft 2 , need not be small quantities. The results 
are equivalent to the Newtonian results for small scale irregularities and p « p c 2 . 

For |i/»| «— c 2 , the Schwarzschild line element is 


With this as an analogy McVittie 53 assumes that the line element is orthogonal 
and of the form (119) for a perturbed Einstein universe, and considers a universe 
of discrete condensations. By supposing in a linear treatment that 4> depends 



■f) 


(dr 2 + r 2 d G 2 ) . 


(131) 
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upon the entire mass of the condensations it was found that the volume of a static 
universe depends on the number of condensations. Modifications 54 were sub- 
sequently made in the treatment which showed that the volume was independent 
of disturbances. 

Bonnor 56 has considered the growth of single condensations in the linear 
approximation for a pressure free fluid and has concluded that in an expanding 
universe the nebulae cannot have originated from infinitesimally small disturbances. 

The most general and elegant contribution to the subject has been made by 
Lifshitz. 46 His basic equations are similar to (110)-(113). In general, coordinate 
transformations are possible which allow four of the h kj to be zero. Which of 
the h k . are made zero determines to some extent the simplicity of the equations 
for a given physical problem. Lifshitz makes the choice h Q . = 0. This choice 
is appropriate for complex disturbances containing rotational motion and gravi- 
tational waves, but the simplest of all disturbances — irrotational motions — are 
concealed within cumbersome equations that require tensorial spherical harmonics 
for their analysis. Lifshitz finds that all perturbations either decay or grow 
very slowly in an expanding universe. 

5.7 Normal Modes of Vibration in Curved Space 

Using r, 0, cp coordinates in 


V 2 0 + k 2 0=0, 


and separating the variables: 0 = '/'(t) 'P(r) V" (&, <fi), we find that only the radial 
function ^(r) is different in the three cases k = 0, ±1. 
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For spherical harmonics of degreen the radial equation is 


,3/2 


(q ,/J r 2 £t) * 

k 2 n(n+l) 

\ dr / 

qr 2 


'P = 0, 


(132) 


where q = ^1 + -i- k r 2 j . We consider briefly the wavefunctions and eigenvalues 


in a space of negative and positive curvature. 
For k = l let 


46, 57 


sin a = 



hence , 



(133) 


and (132) becomes 


1 d 

f . 2 d'P'' 

1 S11T 6 Cl 

!♦ 

j,2 n(n+l) 

_ • 2 « d CL 

V da J 


sin a 

\ / 


L sin z a J 


With <P = n sin 1/2 a as the new variable, this equation is 


1 _d 
sin a da 



K(K + 1 ) - 



n = 0, 


(134) 


and X (K + 1) = k 2 + 3/4. 
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For k = - 1 let 


sinh a = 



hence, da 


dr 



(135) 


and 4* = fl sinh' 1/2 a; then (132) becomes 


sinh a da 


/ . , dn 

( sinh a 

\ da 


k (k+1) - 


( n * i) ! 


sin 


h 2 


n = o, 


(136) 


and k (\ + 1) = k 2 - 3/4. 

Let £ = k 1/2 a, then both (134) and (136) can be expressed in the one equation 


1 _d_ 

sin £ d£ 



\ (k + 1) - 



(137) 


where now k(k + 1) = k k 2 + 3/4, or 


k 


1 . 2 



(1 + k k 2 ) 


1/2 


(138) 


The solutions of (137) are the associated Legendre functions P£(cos £), 
Q£(cos £), and /x = — + n, v - k. 
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Positive curvature (k = 1). In tMs case it is more convenient to set 


k 2 = y (7 + 1). (139) 

Hence, \ +K; and since P£ = P£ , , or P£ - pM , we use v - k . Because 

2 " 'S * 2 

M ± v is an integer, but \x is not an integer, we can use P£ and P'^ as linearly 
independent solutions. 58 These are: 


IT/ n 

y 


- ( tt/2 sin a ) 1/2 P 


2 

2 


+ 


+ 


n 

(cos a) , 
y 


'P; n = (n / 2 sin a) 1/2 P 



(cos a) . 


From the definition 


(140) 



it is seen that as a - 0 (or r 0) 


W ±n - 
y 
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and as n is a positive integer, only is regular at the origin. Therefore 




a <P' n 

n y 7 


a ny ( 7T/2 sin rt 


y/2 p ^ 2 (cos a) , 


i +y 


(142) 


We have assumed that yis integral; it can be shown that this is necessary 
in order that 'I' is periodic or single valued: 

(- cos a) = cos (y-n)7T 'P^, n (cos a). (143) 

Thus the wavefunction is symmetric (antisymmetric) about a = -y w, or r = 2, 
when y - n is an even (odd) integer. This has interesting consequences for ellipti- 
cal space: 0 < r < 2, and spherical space: 0 < r < oo. The transformation: r - 4/ r , 
leaves the metric (73) unchanged.* It is therefore said 59 that elliptical and 
spherical space are indistinguishable because 2 < r < oo is merely a remapping 
of elliptical space. In a perturbed universe, however, elliptical space is not a 
mirror image of 2 < r < 00 for the antisymmetric wavefunctions, and therefore 
it would seem to be an inadequate model of the universe. 

For n = 0, 1, ... y, 


- n 


y M 


s in n a d n + 1 [cos (1 + y) a. 


d (cos a) n + 1 


M n = (1 + y ) 2 [(1 + y ) 2 - 1 ] • • • [(1 + y ) 2 - n 2 ] 


(144) 


But with reversed parity. 
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It follows, for n = 0, 


,po _ sin (1 + 7 ) a 
y (l+y) sin a 


(145) 


and for n = 7, 


- (sin a) y / (1.3.5 • • • • 2y + 1). 


(146) 


Also 


ro-7 

7 + 1 


cos a 


(147) 


U|-7 

7+2 



27+4 

27+3 


sin‘ 



m-7 

7 


(148) 


Hence, the radial functions of the lowest modes are 


V® = cos a , 


¥ 


-1 

1 


1 

— sin a, 
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- sin a cos a, 


V" 2 = — sin 2 
2 15 


a, 


^3 = cos a (1 - 2 sin 2 a), 


1 = 1 sin a ^1 - sin 2 a^, 

W” 2 = — cos a sin 2 a, 

3 15 

V" 3 = _L_ sin 3 a. 

3 105 

In a space of uniform positive curvature the eigenvalues are 

k 2 = y (y + 2), y = 1, 2, 3, • • - 

and m < n < y. The fundamental mode has an eigenvalue of k 2 = X 



Negative curvature (k = -1). For the eigenvalues in this case we use 


k 2 = 7 2 + 1, (151) 

and from (138) X. 12 - - — ± P~ M and are linearly dependent because 
1 2 

M - ~2 +n is half integral, and for the linearly independent solutions P£ and 

P are again chosen. Since F* = pp , we chose v = k,, and therefore 

K 1 K 2 1 

V " = (77 / 2 sinh a) 1/2 P 1/2 + n (cosh a), 
y ’ -i/2+ir 


= ( 7T /2 s inh a) 1/2 P j^ 2 + n (cosh a). 


(152) 


The hypergeometric expression 


p± (1/2 + n) 
- 1/2 + i 7 


(cosh a) = 




- sinh 2 1 
2 



is real, and as a- 0, 


n 

y 


( 1 \ ± (1/2 + n] 

(coth — a) 


(s inh a) 1/2 F ( 1 + 


(3 * -)) 
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and therefore only?"" is regular at the origin. Thus, 


<P” n = a ( 77 / 2 sinh a) 1/2 P 1/2 n (cosh a), 
y n 1 -y v • * > -1/2 + i y 


and because space is open, there are no periodic conditions to satisfy, and 
can have any real value. 

For integral values of n , 


sinh n a dn + 1 ( cos y cl ) 
X d (cosh a) n + 1 


N n = - y 2 {y 2 + 1) (7 2 + n 2 ), 


and for n = 0, 


= sin y a/y sinh a, 


n = 1, 


»P" 1 = (y 2 + 1) _1 (y cot y a - coth a) 


and so forth. 


(153) 


7 


(154) 


(155) 


(156) 
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We have found that the eigenvalue spectra for k = 0, 1, and -1 are 
/<=0: k 2 = y 2 , y 2 >_ 0 , 

k - 1 : k 2 = y (y + 2) , y = 1, 2, 3 • • • , 

k - -1 : k 2 = y 2 + 1 , y 2 > 0 , 

and the lowest eigenvalues are k 2 = 0, 3, and 1, respectively. The eigenvalues 
form continuous spectra for zero and negative curvatures, and a discrete spec- 
trum for positive curvature. 

5.8 Time Dependence of Modes in Einstein’s Static Universe 
From Section 5.1 we have 

(8 ttG p + A) R 2 = 3 (R 2 + k c 2 ), 

(A - 8ttGp/c 2 ) R 2 = 2 RR + R 2 + kc 2 . (157) 

In Einstein’s static model R = R = 0, K = 1, and therefore 

8ttGp + A = 3c 2 /R 2 , 

A - 8ttGp/c 2 = c 2 /R 2 . (158) 
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With the equation of state p= (v - 1) pc 2 , it is found that 


A = c 2 (3^ - 2 )/^Rq • 


( 159 ) 


Equations (125), (123), and (122), become 


0 + c 2 {( v - 1) [k 2 + 2(3^ - 2) v _1 ] - (3v - 2)} Rq 2 0 = 0, (160) 


4ttGR 2 0 Sp = - [k 2 - 3 + 2(3v - 2) v" 1 ] 0 , (161) 


0 = 477G v p cp. 


(162) 


The case of v - 1 has been considered in Section 4.2; we now consider a 
supernuclear model of v= 4/3. Equations (160)-(162) are therefore 


H* 


I 7(7 + 2 ) - 1 0 = 0 , 


(163) 


67tGAc 2 Sp = - 7(7 + 2) 0 , 


(164) 


0 = 167rGpcp/3, 


( 165 ) 


in which k 2 = y(7+ 2) is used. The solutions are 
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Thus all modes of y > 1 are oscillatory, and the fundamental mode y - 1 has a 
solution which grows linearly with time: 

ip = At + Bt -1 . (167) 


If the Einstein universe were static for an indefinite period of time the linear 
growth of the f undam ental mode would count as a serious instability. It is well 
known, 60 however, that the Einstein universe is unstable against perturbations 
in R. Perturbing (157) and using the equilibrium equations (158) we find, to afirst 
order, SR = v\ SR. Therefore, 

SR oc exp (3v - 2) 1/2 ct/R 0 , (168) 

and for v > 2/3, SR grows exponentially. Hence the growth of the fundamental 
mode (167) must be re-examined in a nonstatic universe. 

It was thought 53 that condensations in a static universe would increase the 
volume to V Q + Sv, thus launching the unstable universe on a career of expan- 
sion. Using the line element (119), we find 


SV = (3R 2 /c 2 ) s in 2 a s in 0\p(a, 6, 4>)daddd4>, 


(169) 



in the linear approximation, and therefore §V = 0. Thus, in the linear theory 
condensations do not guarantee that the Einstein universe evolves into an expand- 
ing Lemaitre model. 

5.9 Time Dependence of Modes in Nonstatic Models 

The time dependence of the modes for p « pc 2 has been considered in Sec- 
tion 4, and here we shall consider only the supernuclear model of v = 4/3 and 
A = 0. The equations for pi andS/o are 


i/r' + 5RR -1 <A + (c 2 


/R2 > (§ k2 - 


4xr) 0=0, 


(170) 


477Gc 2 h p - - 3RR -1 ip - R -2 (k 2 c 2 - 3k c 2 + 3R 2 ) 


(171) 


from (128) and (123). With the transformation (93), k = 0, these equations are 


d 2 i// 4 dip 

d x 2 X d X 



4 > - o, 


(172) 


1 8 p 

2 ~ 




' P- 


(173) 


Because y is generally a small quantity (y < 10 "«), these equations are valid 
for k = 0, ±1. In the supernuclear state there is no simple relation between \p 
and hp. The solution of (172) is 
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0 c cy- 3/2 J ±3/2 (sy), 


► 


s = 


(I 


k 2 


- 4^ 


1/2 


(174) 


For the low order modes sy « 1, and for k = 0, ±1 


0 = A + BR" 3 , 


(175) 


where y ot R cc t 1/2 , and 


hp/p - AR 2 + BR -1 , 


(176) 


and A and B are different in the two equations. Thus <b is either constant or 
diminishes in an expanding universe, and 8 p/p increases linearly with time, or 


8 p/p cc R 2 a p“ 1/2 . 


(177) 


In general, 



0 

1 

A / S * n - cos; s 

+ B 

/ * COS S Y VI 

sin Sy + _ i 




\ s x j 

( 

\ s * /J 


(178) 


We consider first flat space: k — 0, s =k//3, and assume that s y >> 1. It follows 
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(179) 


i -2 iky / 1/3 
\p OZy z e A , 

o ^ / iky/ v/3 

op/p oc e A , 

and &p / p oscillates with constant amplitude during expansion and contraction. 

These equations, because of (94), hold for a wavelength ft = Rk " 1 of ft << R 3 / 2 /Rq /2 , 
where R < R n . For a photon or neutrino universe the condition >>1 is simply 

ft << *R - c/(pG ) 1/2 . 

In curved space: k = ± 1 , and for sy >> 1 , 


0o cy 2 exp i 



(180) 


and 8p/p oc0y 2 . These results are the same as (179) in a universe in which the 
supernuclear state holds for y << 1. But in a photon or neutrino universe of 
positive curvature the modes 


k 2 = y(y + 2 ) < 12 , 

or the first two modes 7=1 and 2, grow exponentially. However, since 0 < y <rr , 
their growth is limited. 

We can summarize by saying that the contrast density oscillates with con- 
stant amplitude for short wavelengths and increases linearly with time for long 
wavelengths. 
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6. DISCUSSION 


The vibrations of the universe are moderately well behaved and show no 
signs of a catastrophic growth in amplitude. This conclusion has been reached 
by several authors in various ways. Within the limitations of simple, linearized 
* theory the expanding universe is reasonably stable. This is illustrated in an 

approximate manner by Jeans’ theory. The gravitational frequency is of the 
order i (pG) 1/2 , and the age of the universe is of the order (pG) _1/2 , and hence 
perturbations cannot grow significantly in the time available. 

In the supernuclear state the contrast density oscillates at constant ampli- 
tude for short wavelengths and increases for long wavelengths according to 

Bp/p oc p -1/2 . (181) 

In the subnuclear state the contrast density depends on curvature and the ratio of 
specific heats. The maximum growth rate occurs in a cold universe and 

5p/pap" 1/3 , (182) 

for expansion at zero curvature. The significance of these results depends on the 
magnitude of the initial disturbance that can be legitimately assumed. For arbi- 
trarily small initial disturbances the results (181) and (182) tend to be meaning- 
less. If the disturbances originate from thermal fluctuations, the mean square 
fluctuation in N particles is 40 

((AN) 2 )= kT(BN/Bp) TV (183) 
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at constant temperature T and volume V, where /x Is the chemical potential. The 
right-hand side of (183) is equal to N for a Maxwellian distribution, and is of the 
same order for a relativistic gas of bosons and fermions in which the Fermi 
energy is equal to kT . Therefore, 

S p/p ^ N" 1/2 , 

and this is exceedingly small for large masses containing ordinary gas particles. 

A redeeming feature of the present theory is that the universe in the large 
tends to remain homogeneous and isotropic. If the growth rates (181) and (183) 
adequately accounted for irregularity, then the fundamental modes would also 
increase in the same manner, and in all probability the universe would develop 
pronounced macroscopic anisotropy. What is required is that the universe is 
unstable for an intermediate range of wavelengths which grow rapidly in an 
exponential fashion. At some stage in the expansion the density in various regions 
ceases to diminish with time, and thereafter condensations occur. But such a 
concept demands that the gravitational potential of the disturbance increases 
with time and attains a value of ip ~ GM/A. where M is an island mass and \ its 
radius. But nowhere in the linear analysis is 0 an increasing function of time 
in an expanding universe. 

A linearized theory limited to irrotational motions and gravitational inter- 
actions is open to several criticisms. The neglect of all forms of rotation is a 
gross simplification, since angular momentum is a common and indispensible 
feature of galactic and stellar systems. It seems plausible that at subnuclear 
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density, at least, a treatment based on rotational motions will lead to even slower 
rates of growth owing to the presence of inertial forces. Gamov 61 has proposed 
a primordial turbulent state of large amplitude fluctuations for the initial con- 
ditions. Similarly, Weizhcker 62 assumes an initial state of turbulent gas clouds. 
Bonnor 56 points out, however, that turbulence is more likely to be the result 
rather than the cause of condensations. Furthermore, initial conditions of this 
nature add to the mystery rather than clarifying it, and their postulation falls 
within the province of hypothesis A and renders invalid a linear treatment. 

The assumption that initial disturbances are small, as in hypothesis B, 
demands that rapid growth is possible at some stage in the expansion of the 
universe. This is impossible unless we abandon the rudimentary fluid prescrip- 
tion. Following a speculative vein, we might suppose that because of complex 
fluid properties (such as energy transfer by radiation or neutrinos) the equation 
of state for perturbations is Sp= -(1 - v) §pc 2 and v < 1. In effect, the velocity 
of sound is imaginary. Then, for ^ >>R/R , 

\p oc exp ~ [(1 - v ) c 2 k 2 /GpR 2 ] 1/2 . 

If such an approach to the problem were possible and could yield rapid and large 
growth, then as | \p | — . c 2 , the linear theory breaks down. It is interesting to notice 
that such growth not only could account for differentiation but might also lead to 
the development of phenomena such as the partial concealment of dense masses 
encapsulated in the metric. 63 Such bodies, consisting of matter in a primordial 
state, would then burst forth at subsequent epochs as quasi stellar objects, 64 and 
possibly as explosive nuclei in galaxies. 65 To an observer in such a body the 
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universe is still quite young. These concepts have much in common with 
Ambartsumian's hypothesis , 66 and provide a physical mechanism as its basis. 

Clearly, the origin of structure offers a challenging problem, and from this 
point of view a study of the modes of vibration of the universe has merely em- 
phasized the dilemma that confronts us. If there is no way in which we can 
recover our faith in the efficacy of hypothesis B, then we must turn to hypothesis 
A in the hope that the primordial universe contains structure as the natural con- 
sequence of the properties of matter at very high density. 
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